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Abstract 



We propose a new formula for the adiabatic Berry phase which is based on phase-space 
formulation of quantum mechanics. This approach sheds a new light into the correspon- 
G i dence between classical and quantum adiabatic phases — both phases are related with the 

averaging procedure: Hannay's angle with averaging over the classical torus and Berry's 
phase with averaging over the entire classical phase space with respect to the correspond- 
ing Wigner function. Generalizations to the non-abclian Wilczek-Zee case and mixed 
states are also included. 

>. 
\G , Geometric Berry phase [1] and its classical analog, so called Hannay angle, [2] (see also [3]) 

>-S I have found numerous applications in various branches of physics (see e.g. [4] and [5]). Re- 

\^ • cently, it turned out that adiabatic Berry phase plays important role in quantum computation 

^D . algorithms as a model of a quantum gate in a quantum computer (see e.g. [6, 7]). In this 

paper we present a new formula for the Berry phase which is based on the phase space for- 
mulation of quantum mechanics. This approach sheds a new light into the correspondence 
between classical and quantum adiabatic phases. 

Both Berry's phase and Hannay's angles have been introduced in the context of adiabatic 
!^ '_ evolution in quantum and classical mechanics, respectively. Let us consider for simplicity a 

classical system with one degree of freedom and let the corresponding phase space be param- 
^. eterized by canonical coordinates {q,p). Suppose, that a Hamiltonian H{q,p;X.) depends on 

a set of some external parameters X from the parameter space M and that X are changed 
adiabatically along a circuit C and come back to their initial values, i.e. X(T) = X(0) for 
some T > 0. Now, the classical adiabatic theorem [8] states that the system will evolve on 
5^ I the torus defined by the constant value of the action variable I and the angle variable varies 

according to 

0{T)= [ uj{I;X{t))dt + A0{I;C) , (1) 

Jo 

where the frequency uj{I,X) = dH{I;X)/dI and the additional shift — Hannay angle A9 

— is given by the following integral over an arbitrary two-dimensional region S in A^ such 

that C = as 

AeiI;C) = -^^ II^^^^F%I;X), (2) 

where F^{I; X) denotes the following two- form on A4: 

F=(/;X) = (dxp(/;X)Adxg(/;X)) , (3) 



and (/(/)) denotes the average of /(/, 6) over the torus /. 

Now, let us consider the quantized system defined by H(X.). Clearly, the quantization 
H ^ H is not unique and depends on the ordering of q and p. In what follows we assume 
the Wigner-Weyl (or symmetric) ordering, i.e. for example qp -^ {q.P + PQ)/'^- The quantum 
adiabatic theorem [9] states that a system originally in an eigenstate |n; X(0) ) will remain in 
the same eigenstate |n;X(t) ) with energy £'„(X(t)). Now if the initial state |V'(0) ) belongs 
to the nth eigenspace, then after the circuit C the final state \'4'{T) ) is given by 



|V(r) ) = exp(i7„(C)) exp ( -'- ^ En{Mt))dt ) IV'(O) ) , (4) 



h 



where the Berry phase reads: 

ln{C) = - [ [ F„^(X) , (5) 

and Fn (Berry's curvature) denotes the following two-form on ^A^. 

F^(X) =Im(dxn;X| A|dxn;X) . (6) 

Using semiclassical analysis Berry shown [3] that 

A.(;;C) = -ft5^ = -^, (7) 

where n is considered as a continuous variable according to Bohr~Sommerfeld quantization 
rule / = h{n + ^), with ^ being the Maslov index. 

Both two- forms F^{I; X) and Fni^) live in the parameter space M. Clearly, they are de- 
fined by very different objects: 'classical form' F^{I] X) uses phase-space quantities q{I, 9; X) 
and p{I,9;J^) whereas 'quantum form' Fn(X) is defined in terms of Hilbert space eigenvec- 
tors |n;X). There is, however, equivalent formulation of quantum mechanics which uses 
objects defined on the classical phase space only [10]. There is a direct relation — well 
known Wigner-Weyl correspondence — between functions F = F{q,p) on the classical phase 
space and self-adjoint operators F in the system Hilbert space. If |V') is the state vector, 
then 

{i;\F\i;)= [w^iq,p)F{q,p)dqdp, (8) 



where W^ = W^{q,p) is a Wigner function corresponding to \ip). Moreover, this formulation 
is perfectly suited to semiclassical analysis. It is well known that Wigner function corre- 
sponding to the eigenstate \n) of the Hamiltonian H depends only on / and not on 6, i.e. 
^niQ,p) = Wnil)- In the classical limit, i.e. h ^ 0,n ^ oo such that nh is constant and 
equals Iq, the Wigner function Wn is concentrated on the torus /q 

Wnil) -^ ^6{I-Io). (9) 

If the quantum system depends upon external parameters X which evolve adiabatically then 
Wn is adiabatically constant, or, using the language of the classical adiabatic theorem, Wn 



defines an adiabatic invariant. Now, since Berry's curvature Fn is a measurable quantity it 
may be expressed according to 

F^ = I Wn{I) 'classical quantity' d/ . (10) 

Clearly, this 'classical quantity' has to be related with the 'classical two- form' F'^. Moreover, 
in the classical limit F^ and its classical counterpart F'^ have to be related according to (7), 
that is, 

F„^(X) = -iF'=(/;X). (11) 

The natural choice is therefore 



F^{X) = -27r J Wn{I)F%F,X)dI . (12) 



Example. As an example let us consider a generalized harmonic oscillator [3] defined by 

H{q,p;:K) = ^{Xq^ + 2Yqp+Zp^) , (13) 

where the parameters X = {X, Y, Z) G M^ satisfy XZ > Y"^ (this condition implies that 
the above system describes oscillatory motion round elliptic contours in the two-dimensional 
phase space M^). One shows [3] 

F^U; X) = — ^ iXdY AdZ + YdZ AdX + ZdX A dY) , (14) 

4a; •^ 

with the frequency of the quasi-periodic motion uj = ^J XZ — Y"^. The quantized system 
(according to the Wigner-Weyl correspondence) is given by 

^(X) = i {Xf + Y{qp + pq) + Zp") . (15) 

The eigen-equation Hipn = E^ipn is solved by the following normalized eigenfunctions: 

-iYq^^ 



ipn{q; X.) = y/axn (aq) exp 



2Zh 



where a = y^uj/Zh, and Xn{0 = Nn&~ Hn{C)} with Hn being the nth Hermite polynomial 
and the normalization constant Nn = (2"n!v^)~^'^. Energy levels are given by the standard 
formula En = hio{n + 1/2). The corresponding Wigner function Wn{q,p', X) reads as follows 

Wn{q,p;X) = ^ / dsV'„(g + s;X)Vn(g-s;X)e2W^ 



oo 

oo 



— / dsxn {a{q + s)) Xn {a{q - s)) e2.[(p+y./^).]/ft 



^osc 



WZ{q,p) , (16) 



where W°^'^ is the nth Wigner function corresponding to the standard harmonic oscillator 

with p = p + Yq/Z. Clearly, the X-dependent canonical transformation {q,p) — > (QiP) 
transforms (13) into (17). Now, W°^''^ is given by the well known formula (see e.g. [10]) 

Wn{I) = Wr{q,p) = ^r e~''/^^n(4//^) , (18) 

irn 

where / = H°^^/uj is the action variable and L„ denotes the nth Laguerre polynomial. Finally, 
using 

n + - 

Wn{I)IdI=^^ , 

one finds the following formula for the Berry curvature 

77,-1-1 

F,?(X) = ^F'=(/;X) 

n + - 
= — ^ {XdY AdZ + YdZ AdX + ZX A dY) , (19) 

in perfect agreement with [3]. 

Our basic formula (12) may be generalized in two evident ways: if the classical integrable 
system has A'^ degrees of freedom then the corresponding Berry curvature reads: 

F^{y.) = -{2^f j...jWn{l)F^{l;y^)dl , (20) 

with 

N 
F^(I;X) = ^(dxPfc(I;X)Adx(?fc(I;X)) , (21) 

k=l 

where now one averages over iV-dimensional torus I = (/i, . . . , In)- The second generalization 
corresponds to the non-abelian case developed by Wilczek and Zee [11]. Suppose that nth 
eigenvalue is N times degenerate and let |n, a; X ) (a = 1, . . . , N) span the corresponding A^- 
dimensional eigenspace. Then the Wilczek-Zee curvature is given by the following formula 

fW|(X) = -27r j Wr..,abiI)F%I; X) dl , (22) 

with Wn-ab being the following 'Wigner matrix' 

1 r°° 

H^„;,b((7,p;X) = — / ds(n,a;X|g + s)(g-s|n,6;X)e2^W^ (23) 

Clearly, Wn-^ab is hermitian and hence iF^ G u{N). Now, changing |n, a;X) to |n7a;X) = 
^j,t/a6(X)|n,6;X), with U{X) G U{N), one finds 

fWZ(x) = C/(X)FWZ(x)t/t(X) , (24) 



that is, tensorial rule for the gauge transformation of F^'^. Finally, the formula (12) suggests 
the following generalization for the adiabatic evolution of mixed states. Suppose that p is a 
density operator such that the corresponding Wigner function Wp = Wp{I) is adiabatically 
constant. Following Sjoqvist et. al. [12] one defines a phase of p{T) with respect to p as 
(f) = argTr[f7(t)/j]. Now, (p may be recovered from space-phase quantities as follows: 

cP= f f Fp{X) , (25) 

J JdY:=C 

where the two- form Fp(X.) on the parameter space Ai is defined by 

Fp(X) = -27r / Wp{I)F%F, X) dl . (26) 

Clearly, if the stationary (in the adiabatic limit) state p is pure, then necessarily p = 
|n;X)(n;X| and (26) reduces to (12). 
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